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We present a self-consistent approach to deal with the pairing-fluctuation effects in quasi-two- 
dimensional superconducting systems. Besides the Cooper pairs in the Bose-Einstein condensate, 
there are pairs occupying the excited states, which results in the predominant fluctuations. The 
low-lying excited states are the collective modes. On the basis of ladder-diagram approximation, we 
treat the single particles and the pairs on an equal-footing manner. The Green's function of single 
particles is obtained as an analytic solution to a cubic equation. The bosonic degrees of freedom are 
relevant to the pseudogap physics in the high-Tc cuprates. The superconducting order parameter 
and the transition temperature are substantially reduced from the values of the mean-field theory. 
The calculated phase boundary of superconducting state can reasonably describe the experiment 
data for cuprates. 

PACS numbers; 74.20.-z, 74.62.-c, 74.20. Mn 



I. INTRODUCTION 

The high-temperature superconductors are typical 
quasi-two-dimensional (quasi-2D) electron systems. The 
motion of the electrons is mainly confined in the Copper- 
Oxygen planes with weak coupling between them. In the 
quasi-2D systems, the fluctuation effects may take impor- 
tant role for describing the superconductivity. As stated 
by the Mermin-Wagner-Hohenberg (MWH) theorem,^ 
there is no superconductivity at finite temperature in 
systems of dimensions less than or equal to 2. This is 
because that pairing fluctuations prohibit the electrons 
from coherent pairing. It is therefore easy to understand 
why the fluctuation effects can be significant in quasi-2D 
systems. 

There are a number of approaches dealing with the 
fluctuation effects.^ For the under doped high-Tc cuprates 
(HTC), Emery and Kivelson have argued that the long- 
range classical phase fluctuation of the order parame- 
ter can significantly suppress the transition temperature 
Tc-^ Above Tc, pairing becomes locally without long- 
range phase coherence, especially at the strong-coupling 
regime. On the other hand, the superconductivity can 
be viewed as a consequence of Bose-Einstein condensa- 
tion (BEC) of the preformed pairs. ''^^ Along with this 
approach, much effort has been devoted to investigation 
of the crossover from the weak-coupling BCS supercon- 
ductivity to the BEC of bound pairs. Most of the 
works have been performed for the s-wave pairing be- 
cause of its computational simplicity. 

Here, we note that according to the general theory 
by Goldstone, Salam, and Weinberg, there exist Gold- 
stone modes in the broken-symmetry state of systems 
without long-range interactions. Since the Goldstone 
modes are the lowest excited states for the pairs, the 
excitations of pairs to these states are the most predom- 
inant fluctuations in the system. Before the Goldstone- 
Salam- Weinberg theory, Anderson had studied the collec- 



tive modes in the superconducting state with perturba- 
tion treatment. These collective modes are actually the 
Goldstone modes. The signiflcance of collective modes 
has been observed in the physical properties of other 
systems. An example is the spin- wave theory for the 
Hubbard model at half filling. ^^^^■^ The collective modes 
(spin waves) not only correct the ground-state energy of 
mean- field theory (MET), but also modify the order pa- 
rameter significantly. For the one-dimensional Hubbard 
model, the ground-state energy by the spine- wave theory 
is very close to the exact result; the typical error is about 
1%, but the error of the MET is about 28% at strong- 
coupling limit. Eor the 2D system at the ground state, 
the magnitude of the order parameter at strong coupling 
is reduced to about 60% of the MET value. ^^^^^ Along 
with the perturbation approach, a number of works have 
investigated the pairing fluctuation effects. 

However, for systems of dimensions < 2 at finite tem- 
perature, such a perturbation treatment is not valid as it 
contradicts with the MWH theorem. Even for quasi-2D 
systems at finite temperature, especially at strong cou- 
pling regime, the perturbation treatment is not accurate. 
The reason is that the fiuctuations and the mean-field 
ordering are equally significant. Because the Goldstone 
modes coexist with the coherent pairing and may not be 
regarded as perturbation in quasi-2D systems, an equal- 
footing treatment of the single particles and the collective 
modes is therefore desirable. 

In this paper, we investigate the superconductivity in 
the tight-binding model with d-wave attraction. The 
pairing fluctuations above the mean field stem from the 
excited pairs. According to the many-particle physics, 
we need to find out the bosonic Green's function for the 
pairs. The bosonic Green's function is generally cou- 
pled with that of the single particles. We present a for- 
malism treating the single particles and the pairs on the 
equal-footing manner. In the present approximation, the 
Green's function for the pairs is described by the symmet- 
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rical ladder diagrams. All the Green's functions are self- 
consistcntly determined by a number of coupled integral 
equations. In the BEG picture of the superconductivity, 
the BEG from single pairs begins to occur at Tc- Below 
Tc, with the condensation taking place, the single pairs 
become moving collectively. Even at the ground state, 
there remains the zero-point motion. These bosonic de- 
grees of freedom are relevant to the pseudogap physics 
in cuprates.^''^'^^ We calculate the phase boundary of su- 
perconductivity and compare the result with experiment 
data of HTC. Some results of this work have been pre- 
sented in Ref. 27 



II. HAMILTONIAN 

For describing the electrons, we consider the Hamilto- 
nian 

-H' = X] ^kcl^cka + ;^ X] '"kk'pH^^^ q)p{k', q) (1) 

ka kk'q 

where cj^^ {cka) is the creation (annihilation) opera- 
tor for electrons with momentum-fc and spin-Q, = 
~2t{coiikx + cosfcj,) — 2tz cos kz — n with ji the chem- 
ical potential, Vkki = -JiVkVk' + vt'^t') '^i^h ry^' = 
coskx ± cos ky, p{k,q) = c_k+g/2iCk+q/2t is the pair- 
ing operator, and N the total number of lattice sites. 
Here, we mainly consider the pairing in d-wavc chan- 
nel. By neglecting the s-wave-channel coupling, the in- 
teraction Vkk' is then simply written as Vkk' = —vqkflk' 
with rik = Tjf. For the s-wavc pairing, one just puts 
in r]k = rjf.. For studying the pairing with on site at- 
traction, one can simply set 77^ = 1. In real space, the 
Hamiltonian given by Eq. (1) resembles the t — J model. 
But in the t — J model, double occupation on the same 
site is prohibited. For taking into account of this con- 
straint, we here adopt the approximation introduced by 
Baskaran, Zou, and Anderson^^ assuming the hopping 
integrals t and are proportional to the hole concentra- 
tion 5, e.g., t = tod with to a constant. For the quasi- two- 
dimensional system, t^/t <C 1 is supposed. Such a model 
has been adopted by a number investigators for studying 
the d-wave superconductivity as well as the pseudogap 
phenomena in cuprates.^^~^^'^° Throughout this paper, 
we use the units in which 7i = fcs = 1. 

III. GREEN'S FUNCTIONS 

To formulate our theory, we need the Green's functions 
for the single particles and the pairs. In Nambu's space, 
the Green's function of the single particles is defined as 
in the textbooks by: 

G{k,r-r') = -{TrCk{T)ci{r')) (2) 



where T-^ is the imaginary time-r ordering operator, (• • •) 
means a statistical average, and cj. = (c|.p C-^j^). For the 
Green's functions of the pairs, we firstly need to define 
the structured pairing operators 

Pi(g) = ^^r;fep(fc,g), (3) 



P2{q) = ^J2'^kpHk,-q). (4) 

The operator Pi (g) annihilates a pair of particles of sym- 
metry as described by rjk and of total momentum q, while 

P2{q) annihilates a pair of holes of the same symmetry 
and momentum q. The Green's functions for the pairs 
are defined as 

n^,{q,T-T') = -{TrP^,{q,T)Pl{q,T')) (5) 

These functions are sometimes called as pair suscep- 
tibilities. By the ladder-diagram approximation, after 
Fourier transformation on the imaginary time, we obtain 
the Dyson's equation for 11^;^: 

n(5, Zrn) = x{q, Zm) - vx{q, ZmMq, Z^) (6) 

where n((7, Z^n) and Zm) are 2x2 matrices, Zm = 
i2rmTT with m an integer is the boson imaginary fre- 
quency, and T the temperature, x is the irreducible sus- 
ceptibility. A formal solution to Eq. (6) is obtained as 

U{q, Z^) = [1 + vxiq, Z^T'xiq, Z^) (7) 

In terms of the Green's functions of single particles, the 
elements of x are given by 

X^iiy{q,Zni) = j;^^VkGi,u{kl, Zn + Z,n)G[,fi{k2, z„X 
kn 

(8) 

where fci_2 = k±q/2, = z(2n-|- l)7rT with n an integer 
is the fermion's imaginary frequency, and the subscripts 
/U and = 1, 2 with 1 = 2 and 2 = 1. 

If the BEG of pairs takes place in the system, the boson 
Green's functions n(g, Zm)^s diverge at g ^ and Zm = 
0. This is a general property of the BEG in the boson 
systems without long-range interactions.^"*^ From Eq. (7), 
the divergence means that 

det|l + ^;x(0,0)| =0, (9) 

where x(O)O) is understood as limx(g, 0). (In case of 

BEG, the zero- momentum operators P;^(0) are macro- 
scopic quantities. According to Bogoliubov, they are 
treated as c-numbers.^^'^^ All the boson Green's function 
are defined q ^ 0.) Here, a number of physical mean- 
ings of Eq. (9) need to be noted. Firstly, equation (9) 
is exactly the requirement for the existence of Goldstone 
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modes in the superconducting state. Since the disper- 
sion relation for the pairs is determined by the pole of 
n(g, O), the energy of the pairs Clq vanishes at q ^ 0. 
Therefore, the Goldstone modes are the lowest excita- 
tion states for the pairs and represent the predominant 
fluctuations. Secondly, at Tc, equation (9) reduces to the 
Thouless criterion for the superconducting transition. 

In the Matsubara-frequency space, the Green's func- 
tion of the single particles is written as 

G{k, Zn) = [Zn - CkT3 - Zn)]'^ (10) 

where Zn = Znl with I the 2x2 matrix, a is the Pauli ma- 
trix, and S(fc,z„) is the self-energy. To express the self- 
energy, firstly, we note that the off-diagonal part comes 
from averaged boson fields of momentum g' = 0. At the 
superconducting state, {p{k, 0)) is a macroscopic quantity 
as compared with any other pair fields else. Therefore, 
the predominant contribution is the static mean field 

Si2(A;, Zn) = ;^ l^t'fefe' ipik', 0)) = A^. (11) 

k' 

For our uniform system, we suppose Ak is real. The 
quantity A;; = Ar]k should be difi^erentiated from that of 
the MFT since the fluctuation effect is under consider- 
ation in the present Green's hmction. Secondly, for the 
diagonal part, we take into account of the pair fluctuation 
terms Oof the interaction. By the ladder-diagram ap- 
proximation, the diagonal part of the self-energy is given 

by 

^fifi{k,Zn) = -—'^v^r]l_g/2Gp,p,{k-q,z„-Zm)IitMt,{q,Z, 

qm 

(12) 

Because of the divergence of Hfjfjiq. Zm) at q ^ and 
Zm = as noted above, our treatment for the diagonal 
part of the self-energy takes into account of the predom- 
inant fluctuation effect. With the special form of the 
self-energy, the Green's function can be written as 

G{k, Zn) = [zn - So + (a + S3)a3 + AkO^XjZ, (13) 

Z={zn- So)' - (6 + S3)' - A^, (14) 

where So, 3 = (Sn ± S22)/2, and the arguments {k,Zn) 
have been omitted. 

It should be emphasized here that Eq. (11) is exactly 
consistent with Eq. (9). Any improper treatment of off- 
diagonal self-energy leads to violation of this consistency. 
To see the consistency, we note that Eq. (11) is essen- 
tially equivalent to the gap equation 

A = -^^%Gi2(fc,z„) =-^^77fcAfc/Z. (15) 

kn kn 

On the other hand, the left-hand side of Eq. (9) reads, 
{l + vxn){l + vx22)-v\'i2 = {l + vx+){l + vx-) (16) 



where x's take their values at (g, Zm) = (0,0), xii = X22, 
and x± = Xii ± Xi2- The factor 1 + vx- can be written 
as 

1+VX- = l+^Y.'^k{Gll{k,Zn)G22{k,Zn)-Gl2{k,Zn)]. 
kn 

(17) 

Substituting the Green's functions given by Eq. (13) into 
Eq. (17), we have 

i + v^_ = i + ^J2d/z = o, (18) 

kn 

which is consistent with Eq. (15) and thereby with Eq. 

We need one more equation to determine the chemical 
potential /x. This equation is for the number density of 
electrons, 

— ^Gn{k,Zn)e'-^ = 1-6 (19) 

kn 

where 77 is an infinitesimal positive number. The equa- 
tions (7), (8), (10)-(12) and (19) form the closed system 
that self-consistently determines the Green's functions. 

It is a tremendous task to numerically solve these equa- 
tions because many multi-dimensional integrals over mo- 
mentum and the summation over Matsubara's frequency 
need to be computed in the iterations. However, since 
n(g, Zm) is strongly peaked with a divergence at g ^ 
and Z„i = 0, the diagonal part of the self-energy can be 
approximately given by^^'^**'^^ 

S^^fc, Zn) « -v\lG-^-^{k, Zn)^ Yl 'n^''(9' 

qm 

(20) 

where ^ ' means the summation over q runs small q, and 
the convergent factor e"''^™'' with ai = 1 and a2 = —1 
has been introduced. This convergent factor comes from 
the fact that the Green's hmction Gfifi{k — g, z„ — Zm) in 
the summation in Eq. (12) is connected with the effective 
interaction u^n^^(g, Zm)- The summations over g and m 
in Eq. (20) give rise to a constant 

= -^J2''''li,,{q,Zm)e"''^-\ (21) 

qm 

The constant F is named as pseudogap parameter since 
at Tc there remains a gap in the density of states (DOS) 
at the Fermi energy. Note that the constant is in- 
dependent on the subscript fi because n^y(g, Z„j) = 
^['p.i.Qj ~Zm)- is essentially a measure of the density 
of the uncondensed pairs. By such an approximation, the 
diagonal part of the self-energy is given by 

^ij.ij,{k, Zn) Tl.Gfijj{k, Zn), (22) 
where Ffe = Trjk. 
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With Eq. (22) . wc can get an explicit expression for the 
Green's function of the single particles in terms of F^; and 
A/j. Firstly, note that two equations from the diagonal 
parts of Eq.(13) form the closed system for determining 
the diagonal parts of G{k, Zn)- 



■i- I. 



Go]/Z, 



Gs = - rlG^yz, 



Z = {zn- TlGof - (6 - rlGsf - Al 



(23) 



(24) 



(25) 



where Go, 3 = (Gii ± G22)/2, and the arguments k and 
Zn have been omitted for briefness. From Eqs. (23) and 
(24), one obtains G0/G3 = Zn/^k- Substituting this re- 
sult into one of the equations, e.g., Eq. (23), we can 
obtain a cubic equation for Gq or G3. Instead of writing 
down such an equation for Go or G3, we here introduce a 
function y{k, z„) by Go = Zn{2 — y)/3Tl so that the cubic 
equation for y{k, Zn) and the expression for the Green's 
function look compact. The equation for y{k,Zn) reads. 



(26) 



2 

liri + 2Ai)/{a k) 



where P = 1 + 3(r 



^DKil - zl) and Q = 1 + 
To match the boundary con- 
dition G(fc, Zn) ^1 Zn at n 00, wc should choose the 
real root to Eq. (26). The explicit form of y{k, Zn) reads 



^Q + y/D+ ^Q-VD, D>0 
2\/Pcos((p/3), D <0 



(27) 



where D = - P^, and 



It is then enough to find out y at various x. In terms of x, 
the related quantities can be written as P = 1 + 3(1 — r)x, 
= 1 + 9(r + l/2)x with r = A^/F^. Note that y can 
be expressed by another analytic function, y = f + P/f, 
where / is given by 

^l/2^1/3 



/ = (Q 



(31) 



which is definitely valid at small positive x as seen from 
Eq. (27). For -co < x < 00, equation (31) determines 
the analytic continuation. To extend the definition of / 
in the entire region (—00,00), we write D = 27ax{x — 
Xi){x — X2) with 



xi^2 = —b/a± \J\P' — ar/a, 



(32) 



where a = (r - 1)^, h = (Sr^ + 20r - l)/8. The three 
branch points of D^^^ are xi,2 and 0, with xi < 0, 
a;2|r<i > and a;2|r>i < a^i- Another quantity is 
XQ = — 2/9(2r + 1) at which Q changes its sign: sgn(Q) 
= sgn(a; — xq). Since xi < xq, ReQ is negative at 
X < Xi. Knowing the branch points of D^/^ and the 
property of Q, the analytic continuation of Eq. (31) is 
a straightforward manipulation. Using the polar coordi- 
nates, / = pcxp(i6), we give the expression for p and 6 
in Tables I and II. The angle (p in Tables I and II is given 
by the same functional form as by Eq. (28). 

It is obvious that by setting = 0, one gets the MET 
results. In contrast to the well-defined single particles by 
the delta-singularity in the MET, the single particles in 
the present case have finite lifetimes, and the dispersion 
relation is not defined. The square and cubic roots char- 
acterize the singularity of the present Green's function. 
This leads to the broadened peaks in the density of states 
(DOS) other than the sharp peaks by the MFT.^"^ 



= arccos(g/\/p3). (28) 

The final expression for G{k, Zn) is 

G{k,Zn) = [^„ + 3Afcai/(l + y)+a^T3](2-y)/3r^. (29) 

For more general purposes, we need the formula for 
the retarded Green's function G(fc, to -\- irj) at real fre- 
quency uj. It can be obtained by the analytic continu- 
ation Zn ^ ijJ -\- irj. Firstly, note that G{k, lo + ivj) has 
following property: 

ReGo(fc, — ti> + i?7) = — ReGo(fc, w -|- i??), 

ReGi^3(fc, — w -|- 177) = ReGi,3(fc, w -I- zry), 

ImGo(fc, — w -|- i?7) = ImGo(fc, a; + iry), 

ImGi,3(fc, —'^ + ill) = — IiiiGi,3(A;, w + ijy), 

Therefore, we need only to perform the analytic continu- 
ation for a; > 0. Secondly, y can be regarded as a function 
of '^l/i^l - Zn)\zn^u,+iri -^x + ir]sgn{uj) = X with 



X = 



p2 
^ k 



— CO 



2 ■ 



(30) 



IV. TRANSITION TEMPERATURE 

In this section, we investigate the pairing-fiuctuation 
effect in the transition temperature Tg. Because the fluc- 
tuations play important role in quasi-2D systems, Tc can 
be substantially suppressed from the values by the MET. 

At Tc, A = 0, from Eqs. (13) and (22), we can imme- 
diately obtain a more explicit expression for the Green's 
function 



G{k,Zn) = iZn+^k(^3)(i 




(33) 



The properties of pairs are described by the suscepti- 
bility. Firstly, we analyze the behavior of the irreducible 

susceptibility (at Tc) x{q, Zm) = Xo(<?, Zm)+X3{q, Zm)cr3 
at small q and Z^- The expressions for the Pauli com- 
ponents are given in the Appendix. Note that xo{q, Zm) 
and xsili Zm) are respectively even and odd functions of 
Zm- At Tc, equation (9) means 1 -|- t'Xo(O)O) = 0. At 
small q and Zm, we have 

1 + vxo{q, Zm) « cig2 + c^q^ + 0{Zl), (34) 
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vx3{q,Zm) ~ -dz„ 



(35) 



where ci, Cz, and d arc constants, q in the right-hand- 
side of Eq. (34) represents the in-plane components of 
momentum. The Qz term comes from the ^;-direction mo- 
tion of the particles. To the leading order, the constant 
Cz is proportional to {tz/tY (see the Appendix) which is 
much less than ci. The susceptibility Jli\{q,Zm) is then 
approximated by 

vnr^{q,Zm) « -l/{cxq^ +Czql - dZ^). (36) 

The pole of IIii (g, Q) gives rise to the low-lying excitation 
energy = (cig^ +Czql)/d for the pairs. The constants 
d/2ci and d/2cz represent the effective masses of in-plane 
and out-plane motions, respectively. 

According to Eq. (21), the pseudogap parameter F is 
calculated via 

=-iEVnn(g,^„)e^'"'' 



qm 



^ n„-2 



Tv 

Nd ^ Cta-Zn 
qm 



(37) 



where B{^q) = l/[exp(f2g/T) — 1] is the Bose function. 
Equation (37) clearly reflects the fact that the pairs are 
independent on each other. So the total number of the 
pairs is the summation of the boson-occupation numbers 
on each state. A character of the independent pairs (sin- 
gle pairs) is fig cx at qz =0. The g-integral in Eq. (37) 
is over a region of small momentum. Since depends 
weakly on the out-plane wave number qz, the integral 
over qz can be taken in the range (— tt, tt). Though the z- 
direction dispersion is weak, it prevents the summation 
over q from a logarithm divergence at the g = limit 
and ensures a finite transition temperature Tc- This is 
in agreement with the MWH theorem. The cutoff qc for 
the in-plane wave number is determined such that the 
largest in-plane energy = 2r since a pair that is a 
bound state of two particles is meaningful only within 
the gap. Note that Q.q is a two-particle excitation en- 
ergy. The largest gap for a pair of particles should be 
4r. But regarding the d-wave pairing, the average gap is 
less than 4r. The cutoff 2r is a reasonable choice. We 
have examined the insensitive dependence of and T on 
this cutoff. Even by putting in qc = it, the results for 
and r change very little. 

Equations (9), (19) and (37) self-consistently deter- 
mine the quantities Tc, /i, and F. The result for Tc as 
a function of hole concentration S (solid line) is plot- 
ted in Fig. 1. In the numerical calculation, we take 
v/2tQ ~ 0.1 and tz/t ~ 0.01 for describing the cuprates.^'' 
For La2CuO, v ~ 0.13 eV has been determined by 
experiments.^^ Therefore, the choice ofv/2to corresponds 
to to ~ 0.65 eV, which is consistent with estimates from 
experiment data."^'' The maximum transition tempera- 
ture Tc,Max ~ 0.01575*0 (« 118 K) obtained by the 



present theory appears at a certain S between 0.125 and 
0.15. The MET result and the experiment data"^^ are 
also shown in Fig. 1 for comparison. The phase bound- 
ary of superconductivity given by the present theory is 
close to the experiment results. In the optimally to over 
doped region, the present theory fits the experiment data 
very well. The improvement of the present theory to the 
MFT is significant. Because of the fluctuations, Tc is 
substantially suppressed from that of the MFT. Espe- 
cially in the under doped region, contrary to the MFT, 
Tc of the present calculation decreases with decreasing 
S and vanishes at 5 = 0. Tc function of 5 by the 
present calculation shows a parabolic behavior the sim- 
ilar as the experiment data, while the MFT result is a 
monotonically decreasing function of 6. 

By the MFT, at T = < Tc,mft, there are pairs in 
the condensate. When pairing fluctuations set in the sys- 
tem, all these pairs are changed to be incoherent pairs. 
This results in the reduction of T^. This fluctuation ef- 
fect is more notable at the strong coupling regime where 
the strength v/toS is large; the particles intend to pair 
independently. This is clearly reflected by the pseudogap 
parameter F shown in Fig. 1 as the dashed line. 

One may ask why Tc^mft 7^ at (5 = 0. This is because 
there exists local pairing at (5 = by the MFT. We can 
analytically show that Tc,mft is finite at 5 = 0. Consider 
the gap equation at Tc 



^^ry^anh(a/2re)/a 



At 5 = 0, we have = 0, and thereby tanh(^fe /2Tc)/^k 
l/2Tc. One then obtain 

T - — Vr,2 _ 1 



For our system, v = 130 meV, Tc^Max w 0.01575*0 = 
10.2 mcV, therefore v/ATc^Max ~ 3.2. At small S, £,k is a 
small quantity, and ta.nh{S^k/2Tc)/S,k can be expanded as 
tanh(^fe/2Tc)/^fc w l/2Tc - ^1/12T^. The equation for 
Tc reads 



24r27V 



Apparently, Tc < w/4 at small 5. 

In passing, we here explain the meaning of pseu- 
dogap parameter using the Green's function given 
by Eq. (33). We consider the spectral func- 
tion at Tc, A{k,E) = -lmGn{k,E + iri)/-K = 



E^)/{E^-ek){E + ik)l2T^Tl, which is 



nonzero only for E"^ — 4r^ < ^ < E"^ . The noninter- 
acting delta-function peak becomes a square root singu- 
larity. Near the Fermi energy, the fc-spacc is constrained 
so that the volume decrease at ^ 0, resulting in a sup- 
pression of DOS at the Fermi energy. This leads to the 
formation of a pseudogap in the DOS.^'' 
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V. COLLECTIVE MODES AND PSEUDOGAP AT 
THE SUPERCONDUCTING STATE 

At the superconducting state, the properties of the 
pairs are qualitatively different from that we have dis- 
cussed in the previous section. The expression for the 
pseudogap parameter will be also different from Eq. (37). 
Below the transition temperature, with the BEC taking 
place, the uncondensed pairs become moving collectively. 

To investigate the properties of the pairs, we start from 
analyzing the susceptibilities. The low-energy pairs of 
particles are described by nii(g', Z^) at small q and Z^- 
This function is given via 



villi = 1 - [1 + v{xo - X3)]/D, 



D = {l+vx-){l + vx+)-v'xi 



(38) 



(39) 



where x± = Xo ± Xi, the arguments (g, Zm) have been 
omitted for briefness. (The denominator D should not 
be confused with the quantity D in Eq. (27).) Since 
1 + vx-{0,0) = as shown by Eq. (18), at g ^ and 
Zm — > 0, we have 



1 + vx-iq, Zm) ~ cig^ + c^g^ 



vX3i<l,Zm) ~ -dZ„ 



cqZ^, 



(40) 



(41) 



The constants ci and Cz are different from that we ob- 
tained in Eq. (34) because of a contribution from xi in 
the present case. Note that 1 + vx+ = 1 + vx- + 2i'Xi , 
we get 



l+VX+iq,Zra)~2p + c[q' + c'^_qi, (42) 
where p = vxi{0, 0). For the denominator D, we obtain 
D = {ciq^ + c,ql){2p + d^q^ + <gf ) - (2pco + rf^)^^ 

= (2pC0 + rf2)(fi2_^2^), 

(43) 

which defines the excitation energy fig of the pairs. At 



besides the constant terms in Eqs. (42) and (44) in order 
to formally reproduce the same formula for flu as at T^. 

Similarly, we can obtain the formulae for ni2, 1121, and 
1122. The form for the matrix 11 then is 



n{q,Zm) = 



M+{q) M-{q) 



(46) 



where M±{q) are two matrices. To the leading order 
0{l/^q), we have 



M+(g)=M_(g) = 



-(1-ai). 



(47) 



Equation (46) reminds us that the pairing operators 
PHq) = [Pl(q),Pl{q)\ can be represented by two eigen- 
modes Ci(g), and C2{q)- 



P{q) = A,Ci{q) + B,Cl{q), 



(48) 



where Aq and Bq are two spinors, determined by AqA'l^ = 
M+{q), and BgSj = M_(g). To the leading order 
0(1/ we get 



^5 = ^J = >/s:si;<^'-''- 



(49) 



These two eigenmodes are the collective modes, with the 
same energy 0<j. 

By applying Eq. (45) to Eq. (21), the pseudogap pa- 
rameter r is now calculated via 



__Tv ^,p + c'^q + dZ, 



Zi, 



(50) 



Carrying out the Matsubara sum, we have 



Tv 



l. P + Cq 
2} fla 



- d/2}, (51) 



where the cutoff qc for t he in-pla ne wave number is now 
determined by Qq^ = 2\/A^"+r^ by the similar reason as 
for Eq. (37). The summation over q mainly comes from 
the first term [B{Qq) + ^]p/Qq because which diverges 



very small q, ^q « with u = 2pco + q-^ at g ^ 0. Clearly, equation (51) is different from 



d . At qz =0, flq (X q a.t q ^ 0. This is one of the 
features of the collective modes, in contrast to that of 
the single pairs flq oc q"^. By the same consideration, one 
can approximate 1 + vxo as 

l + vxo{q,Zm)^P + c°iq'' + clql (44) 

The function vllii{q, Zm) is then approximated by 

lp + c° + dZ„ 



vUii{q,Zm) 



u Zl-ni 



(45) 



where c° = c^q^ + c°q'^. Note that at Tc, p = 0, c? = 
c[ = ci, and = = c^. We keep the small terms 



Eq. (37). The total number of pairs cannot be written as 
the summation over the Bose distributions because the 
pairs move collectively. We note that besides the boson 
occupation number B{ilq), the term 1/2 represents the 
contribution from the zero-point motion. This is another 
feature of the collective modes. At the ground state, 
B{^lq) = 0, the integrand in Eq. (51) behaves like 1/q 
at (/ — > 0. Two conclusions are deduced from this fact: 
(1) is finite at the ground state. (2) Broken-symmetry 
state may survive in a pure 2D system at T = 0. This 
is in agreement with the perturbation theories. The 
numerical results for T have been shown in Ref. 27 

Why can the collective modes appear below Tc? Con- 
sider the case of generating a pair of particles of total 
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momentum q in the system. The generation process in- 
fluences the condensate. This pair comes from a superpo- 
sition of the fluctuation of all the pairs in the condensate. 
This is not the case at T > Tg. There is no condensation 
above Tc. Each pair can be constructed only from two 
single particles. So the pairs with difiierent momentum 
are independent with each other. 



VI. CONCLUSION AND REMARKS 

In summary, we have studied the pairing fluctuation 
effects in the quasi-2D superconducting system. The fluc- 
tuations over the mean-field pairing come predominately 
from the Goldstone modes. We treat these collective 
modes and the single particles on the equal-footing man- 
ner. The respective fermion and boson Green's functions 
for the single particles and the pairs are self-consistently 
determined by a number of integral equations. The 
single-particle Green's function is given explicitly by Eqs. 
(27) and (29). The pairing fluctuations result in lifetimes 
for the single particles. 

We have investigated the superconductivity in cuprates 
using the tight-binding model with d-wave attraction. 
is substantially suppressed from its MFT value. The 
phase boundary of superconducting state given by the 
present theory can reasonably describe the experiment 
results. Also, the pairing fluctuations are relevant with 
the pseudogap physics in cuprates. 

The form of the Green's function given by Eq. (33) 
may be still meaningful at T slightly above Tc- It is 
obtained by observing that the long-range fluctuations at 
Zm = have divergently contribution to the self-energy. 
At T > Tc, though there is no such a divergence, the 
contribution comes from the long-range fluctuations at 
small Zm is still large. The experimental observations 
indicate that the pseudogap parameter weakly depends 
on T slightly above Tc-^^ Therefore, the Green's function 
obtained at may be useful for studying the properties 
of the pseudogap state near Tc. 

For the phase boundary, there is still obvious discrep- 
ancy between the present theory and the experiment at 
very under doped region. This may date from the crude 
treatment of the short-range pair correlations. Local 
pairing without long-range phase coherence is not fully 
taken into account in the present model. Besides this, 
the short-range antiferromagnetic coupling is not cor- 
rectly counted in. To describe the antiferromegnctism 
in cuprates at very small i5, one needs to restart with the 
t — J model. 
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In this Appendix, we give the expressions for the 
Pauli components of the irreducible susceptibility x = 
Xo + Xi^i + Xs'^s- Then, we derive the equation for the 
coefficient describing the ^-direction dispersion rela- 
tion of the pairs. 

The Pauli components of x are given by: 



{q,Zm) = ■j;^^r]k[Goiki, z:^)Go{k2, Zn)-G3{ki, Z+)G3{k2, Zn) 



N 



kn 



{Al) 



Xi{q,Zm) = ^Y.VkGiik2,z+)Gi{ki,Zn), {A2) 

kn 



X3{Q,Zm) = — "^nlGsikl, Zn)[Go{k2, z-)-Go{k2, z+)], 



kn 



(A3) 

where ki^2 = k ± q/2, and = Zn i Z^- Apparently, 
Xo{q,Zjn) and xiil, Z^) are even about 

Zm ^ Zm^ 

while xaiq, Zm) is odd. 

The quantity in Eq. (40) is defined as 



=f4x-(9.,0)U.=o 

= -1^ E vlG^{k,Zn)i,G^{k,Zn)s^, {A4) 

knfx ^ 

where /i = 0, 1, and 3, so = —1, si = S3 = 1. For our 
quasi-2D system, the dependence of the Green's func- 
tions on k^ is only through = — 2t(cosfcj. + cosfcj,) — 
2tz cos fc^ — /i, the derivative of the Green's functions in 
Eq. (A4) can be calculated by 

-^Gij,{k, Zn) = 4:tl sin^ k^—-^G^{k, Zn)+2tz cosk^—G^{k, Zn). 

OK^ Ol;f. Ol;k 

(A5) 

Moreover, t^^/t <C 1, to the first order of tz/t, we have 



d - 

Gu,{k,Zn) = G^{k,Zn) - 2tzCoskz — Gf,{k,Zn), {A6) 

Ot,k 



where the bar functions mean that in which is set to 
0. Substituting (A5) and (A6) into (A4) and earring out 
the fcx-integral, we finally get 

^^ = ^12'jli^G,{k,Zn)rs,. iA7) 

The fc-summation in Eq. (A7) is essential a two- 
dimensional in-plane fc-intcgral. 

By the similar procedures, we can get the equations for 
c'j, appeared in Eq. (42) and c° in Eq. (44). The equation 
for c'z is obtained from Eq. (A7) by setting si = —1. For 
c°, we set si = 0. 
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TABLE I. 
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— pexp{i9) 
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TABLE H. 


Function /(x 4- 177) 


= pexp(i6 


at r > 1. 
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FIG. 1. Transition temperature Tc as a function of hole 
concentration 5. The sohd and dotted Unes represent the re- 
sults of present approach and the MFT, respectively. Both 
theoretical calculations use the same normalization con- 
stant Tc,Max = 118 K. The symbols indicate the experi- 
ment data for cuprates:^® Yi_xCaxBa2Cu306 (solid squares), 
Yo.9Cao.iBa2Cu307-y (open squares), La2-xSrxCu04 (open 
diamonds), Yi_xCaxBa2Cu306.96 (solid diamonds), and 
YBa2Cu307-y (open circles). The result for the pseudogap 
parameter F is shown as the dashed line. 
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